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ABSTRACT 

We derive the fundamental parameters (temperature and luminosity) of 107619 Hip- 
parcos stars and place these stars on a true Hertzsprung-Russell diagram. This is 
achieved by comparing BT-Settl model atmospheres to spectral energy distribu- 
tions (SEDs) created from Hipparcos, Tycho, SDSS, DENIS, 2MASS, MSX, AKARI, 
IRAS and WISE data. We also identify and quantify from these SEDs any infrared 
excesses attributable to circumstellar matter. We compare our results to known types 
of objects, focussing on the giant branch stars. Giant star dust production (as traced 
by infrared excess) is found to start in earnest around 680 L . 

Key words: stars: fundamental parameters — solar neighbourhood — stars: mass- 
loss — circumstellar matter — infrared: stars — Hertzsprung-Russell and colour- 
magnitude diagrams 
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1 INTRODUCTION 

Spectral energy distributions (SEDs) have long been the pri- 
mary method of understanding stars. Colour-magnitude di- 
agrams, which can be quickly made from photometric data, 
enable one to explore various facets of stellar populations, 
such as stellar mass and evolutioanry state. However, these 
does not present the information at its most basic physi- 
cal level: the stellar temperature and luminosity. These rep- 
resent the fundamental ideals of stellar modelling, and are 
theoretically free from biases introduced by photometric cal- 
ibration, interstellar reddening and similar phenomena. 

While transformations to these parameters can be 
achieved through colour-temperature relations and bolo- 
metric corrections, these are limited in scope. Most impor- 
tantly, the wavelength coverage of the observations means 
that well-defined solutions do not always exist for these re- 
lations (e.g. for very red stars, or for observations only cov- 
ering wavelengths longer than the SED peak). Using the 
entire wavelength coverage available allows better determi- 
nation of stellar temperature when a wide temperature range 
is present among a stellar sample. This also allows finer con- 
trol of data quality. All-sky surveys are, in particular, prone 
to contain some poor-quality data due to the large flux range 
they are required to cover, which leads to the saturation of 
bright sources, and the volume of data, which limits the abil- 
ity to match photometric routines to particular situations 
(e.g. in areas of high stellar density or nebular emission). 
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Stellar variability can also cause improper colours to be re- 
ported, which can be reduced by using multiple epochs or, 
equivalently, multi-wavelength data. In this manner, we can 
provide more-robust estimates of parameters for individual 
objects, allowing them to be placed on the true, physical 
Hertzsprung-Russell (H-R) diagram. 

Perhaps the greatest benefit, however, is the ability to 
detect excess flux at a particular wavelength, by providing a 
reference model flux against which fluxes in individual pho- 
tometric filters can be compared. This is particularly helpful 
in the infrared, where colour-magnitude diagrams based on 
only part of the SED can fail to identify sources exhibiting 
emission in addition to the stellar photosphere. Predomi- 
nantly, these sources are either very young stars (pre-main- 
sequence T Tauri stars or Herbig Ae/Be stars), rapid rota- 
tors (classical Be stars), or evolved stars. This latter group 
is mostly comprised of mass-losing red and asymptotic giant 
branch (RGB/AGB) stars, on which we focus our discussion. 

Previously, only colour-magnitude diagrams have been 
used to interpret o ur closest stellar neighbours (e.g. 
iPerrvman et al.l 1 19951 ') . We are no w able to take the data 
returned by the Hi pparcos satellite ()Perrvman fc ESAlll997l : 
van Leeuwepi r2007) and match it with other all-sky surveys 
to produce a true H-R diagram of the local Solar neighbour- 
hood. 

In doing so, we can identify and characterise stars with 
weak infrared excesses which may be otherwise missed by 
convention al colour cuts. While this has bee n attempted 
previously l|lta et alj|20ld : iGroenewegenl [20121 : we later dis- 
cuss these papers in context), this work represents the first 
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by higher-resolution, greater-sensitivity WISE and AKARI 
data; also DENIS photometry was replaced by 2MASS pho- 
tometry, which is less prone to saturation. These substitu- 
tions were performed as follows: 

• Wi and W2 replace MSX B\ and B2, respectively; 

• AKARI [9] replaces MSX A; 

• W 3 and/or IRAS [12] replaces MSX C and D; 

• W 4 and/or AKARI [18] replaces MSX E; 

• 2MASS J and K s replaces DENIS J and K s ; 

• IRAS FSC data replaces IRAS PSC data. 

SDSS bands were dropped under certain conditions, namely: 

• u was dropped if u > Bt + {Bt — Vt) + 2 mag; 

• g was dropped if g > Bt + 2 mag; 

• r was dropped if r > Vt + 2 mag; 

which correspond to ranges beyond which the SDSS data 
(which is prone to saturation for the Hipparcos stars, most 
of which are comparatively bright) cannot be matched in 
simultaneity with the Tycho data to any stellar model. 

The resulting initial input catalogue contains 109 661 
Hipparcos stars with data from it-band to 25 /im. The source 
statistics for this combined catalogue are listed in Table [TJ 



2.2 Fitting the SEDs 

Th e SEDs were fit te d usi ng the code initially described 



in iMcDonald et alj (2009) and modified in the subs e- 
ouent naners jBover et all 2009al: IMcDonald et alj|2010allbl: 



time such a process has been applied to the entire Hipparcos 
dataset and in the context of the stars' absolute, fundamen- 
tal parameters. 



2 FUNDAMENTAL STELLAR PARAMETERS 
2.1 Input data catalogue 



The new Hip parcos (H p ) / Tycho (Bt, Vt) reduction 
( van Leeuwcn 2007) was used as the primary astrometric 
and photometric catalogue, to which the other catalogues 
were matched. Additional data was sourced from the follow- 
ing surveys: 

• Sloan Digital Sky S urvey (SDSS-III) Data Release 8 
(|Aihara fc et alJliullh : ugr-band optical data were in- 
cluded. These data are heavily affected by saturation for 
the Hipparcos sample. The iz-band photometry were left 
out entirely, and bad data from the ugr-bands were iden- 
tified as described in H2.3I 

• DENIS Consortium 2005 data releas^U *Gunn JK B -band 
data were included for sources where 2Gunn > 9-5 mag. 
JA's-band near-IR data were used only when 2MASS data 
were not available (see below). 

• 2MASS All -S ky Catalog of Point Sources 
(|Skrutskie et al.l l2006h : JHK B -ba.nd near-IR data 
were included. All photometry was initially included, 
regardless of data quality: bad data was later removed as 
described in ij2,3l 

• MS X Infrared Astrometric Catalog (Eg an fc Pried 
1996]), incorporating six bands (Bi, B2, A, C, D and E) 
covering 4-18 /im. Note that, at ~18" resolution, these 
data exhibit problems from source blending. 

• The AKAR I-Hipparcos cross-correlated catalogue 
(|lta et al.ll201Ct ). covering 9 and 18 /im, hereafter AKARI 
[9] and [18]. 

• IRAS catalogue of Point Sources, Version 2.0 (PSC; 
iBeichmann et al1ll988|). and Faint Source Catalog, |6| > 
10, Version 2.0 (FSC: lMoshir fc et al.|[l990l l. both limited 
to the 12- and 25-/im bands, hereafter IRAS [12] and [25]. 
At «l' resolution, these data also exhibit problems from 
source blending. 

• Th e WISE Premiminary Data Release l|Wright et al.l 
|2010D . incorporating four bands (Wi through Wi) cover- 
ing 3.5-22 /im. 

As the surveys cover a large range of wavelengths and 
have very differently-sized point-spread functions (PSFs; Ta- 
ble [TJ), one must be conservative in declaring two detec- 
tions as arising from the same source. DENIS, 2MASS and 
AKARI detections were considered to be matched with Hip- 
parcos detections when their co-ordinates matched to within 
1", SDSS and WISE data were matched when within 3", and 
MSX and IRAS data were matched within 5". 

Due to the differing depths and similar wavelength cov- 
erages of the near- and mid-IR catalogues, substitutions 
were made to choose only the most discerning data. Specif- 
ically, where possible, MSX and IRAS data were replaced 

1 VizieR On-line Data Catalog: 11/263 



IWoods et al]|201ll : IMcDonald et al.ll2011bl ldldh. This code, 
hereafter referred to as Getsed, is optimised to detect low- 
contrast infrared excess arising from circumstellar dust. We 
have made some further revisions to the code as detailed 
below. 

Getsed works by minimising the x 2 statistic between 
the observed photometric data and a set of synthetic stel- 
lar spectra to determine stellar temperature and luminosity. 
This requires user-defined stellar mass, metallicity and dis- 
tance, and an interstellar reddening, and given appropriate 
filter transmission curves. A grid of temperatures between 
2400 and 60 000 K is set up, in steps of 400 K. Blackbodies 
of these temperatures are fit to the dereddened photometric 
data and the model with the lowest \ 2 chosen. A finer grid 
is set up and the blackbody temperature is iterated to the 
nearest 25 K. This temperature is then used to derive the 
stellar luminosity and surface gravity to first order. 

For this work, we do not know the mass or metallic- 
ity of our stars, which are required to select the correct set 
of synthetic spectra. Nor do we know the interstellar red- 
dening toward the stars. We assume that the metallicity 
is solar and that the interstellar reddening is zero. An as- 
sumption of solar metallicity is reasonable for nearby Galac- 
tic sta rs: the true values sc atter around the solar value of 
[Z/H] l|Luck fc Heiterl [20051 ) . with the scatter imparting a 
<3% error to the temperature fit for the majority of stars. 
The magnitude of this error is similar to that imparted by 
good-quality photometry. 

A larger error is imparted by interstellar reddening, 
which makes the star appear dimmer and cooler than it 
actually is. This can be significant in the case of distant 
objects, or those in the Galactic Plane. In practice, stars 
which suffer from significant interstellar extinction tend to 
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Table X. Number of sources used from each catalogue. 



Catalogue 


Wavelength 


Beam size 




Sources used 






(jnn) 


(") 


Original 


Combined 


Final 








data 


catalogue 


catalogue 


Hipparcos 


0.528 


~0.5 


117 956 


109 661 


107 616 


Tycho 


0.420, 0.532 


~0.5 


118 924 


109 624 


107 586 


SDSS 


0.354-0.623 


>0.5 


32 253 


30 368 


27 420 


DENIS 


0.786-2.20 


>0.5 


60 083 


2 856 


2 762 


2MASS 


1.25-2.20 


>0.5 


104 324 


104 297 


104111 


MSX 


4.29-21.3 


RJ18.3 


7 663 


7 336 


3153 


AKARI 


10.5, 18.4 


2.4, 2.3 


48 078 


48 013 


47762 


IRAS 


12, 25 


106, 106 


19 728 


16 001 


15 533 


WISE 


3.35-22.1 


5.8-11.8 


64192 


64102 


63 883 



Notes: The ground-based optical beam sizes of SDSS, DENIS and 2MASS are limited by seeing, 
hence presented as approximate lower limits. Beam sizes of scanning satellites are not circular: 
here, the equivalent-sized circular aperture is given instead. The combined catalogue contains all 
Hipparcos objects for which we were able to obtain any matching infrared data. The final cata- 
logue contains only the objects remaining after the data quality cuts described in the Appendix 
had been carried out. 



be the brighter giant stars, which also suffer from significant 
parallax errors. Bright giants are often subject to radial pul- 
sations. These can change the fractional contribution of cool 
and hot spots on the stellar surface to the star's total light, 
leading to changes in the astrometric c entre of light. Thes e 
can impart substantial parallax errors (|van Leeuwenfe 007). 
which can be sufficiently large that extinction is not the pri- 
mary source of error in the placement of these stars on the 
H-R diagram. 

As we do not know the mass of individual stars, we 
estimate it from the best-fit blackbody temperature and lu- 
minosity. As this only affects the stellar gravity, which has 
a minimal effect on the overall SED, we need only to ap- 
proximate the actual mass. We estimate the stellar mass by 
assuming each star is either a main-sequence star or a giant. 
Giants are determined to be stars with: 

( (2.25 x 10~ 4 T) 7 ifT> 6000 K 
\ (6.5 x 1CT 4 T) 7 otherwise 1 

where T is the determined effective temperature in Kelvin 
and L is the determined luminosity in solar units. 

For main-sequence stars, we use a mass-temperature re- 
lation based on a s olar-metallicity, ze ro-age main-sequence 
(ZAMS) isochrone dDotter et alj|2008h . We cannot estimate 
the mass of giant stars so easily. Most stars below the RGB 
tip (L » 2500 Lq) will be the more-numerous, older stars of 
~1 M.0. More massive giants survive to much higher lumi- 
nosities on the AGB, thus we expect very luminous giants to 
be consider ably more massive. Based on the aforementioned 
isochrones |Dotter et alj [20081 ). we assume a mass for the 
giant stars of M = (L/25OOL ) 2/3 M , with limits placed at 
1 and 20 M . 

The first-order determination of temperature, lumi- 
nosity and stellar gravity (from the black-body fit), are 
used as the initial parameters for our synthetic spectra. 
Previously, we have u s ed th e marcs mo del atmospheres 
[ Gustafsson et al.lll975l , l200Sh described in lMcDonald et all 
(2009). For this work, however, we instead use the BT- 
Settl models of lAllard et alj (|2003f ) . as the temperature 
spacing of the grid models is finer in the 4000-6000 K re- 
gion. We have found this to reduce artifacts caused by in- 



terpolation between grid points for stars with poor-quality 
photometry. 

Getsed takes the model spectra grid and performs a 
linear interpolation in temperature and surface gravity (and 
metallicity, if required) . The synthetic spectrum is convolved 
with the observed filter transmissions and reduced to a set 
of expected photometric fluxes. These are then normalised 
to the observed photometric fluxes (the constant of normal- 
isation determining the luminosity) and a x 2 is determined. 

This process is first performed on the temperature grid- 
point immediately cooler than the blackbody temperature. 
Getsed then calculates the x 2 f° r the neighbouring tem- 
perature gridpoints and continues until a x 2 minimum is 
detected. The temperature corresponding to the \ 2 mini- 
mum is used as a new starting point, a new surface gravity 
is calculated, and \ 2 i g determined for 128 K steps between 
the neighbouring models. A new \ 2 minimum is determined, 
the temperature step is halved, and the process re-run until 
the temperature is fit to within 1 K. The calculated stellar 
effective temperature, luminosity, surface gravity and model 
photometric fluxes are written to disk. 

As with our previous uses of this code, we have only fit 
photometric data with short wavelengths (here we require 
A < 8 /im). Circumstellar dust will still cause some opacity 
in the optical, but the obscuration must be relatively small 
for it to have been observed with Hipparcos, and absorption 
of optical flux is relatively easy to identify when it is re- 
radiated in the infrared. 

Stars which are heavily extincted will have SEDs that 
become double-peaked. In these cases, Getsed will not be 
able to fit a mod el spectrum to it. T he most-extincted stars 
(e.g. IRC+10216: lKwan fc Hi"iHll977l ) may be sufficiently op- 
tically obscur ed that they do no t feature in the Hipparcos 
catalogue (cf. iBover et al.ll2009bl ). This becomes important 
in the removal of bad data (see Section I2.3|l and we remind 
the reader that our H-R diagram is therefore incomplete, 
even at high luminosity. 
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2.3 Removing bad data 

Each star now has a series of photometric datapoints for 
which an observed and a modelled flux is known. The ratio 
(R) of observed/modelled flux therefore gives the excess or 
deficit flux in that band. We give this as a pure ratio, rather 
than in terms of an a-based excess, as we do not include the 
errors in the photometric data in our model. While this may 
appear surprising, the reported errors on photometric data 
are almost invariably much lower than the absolute error be- 
tween catalogues. Such 'bad' data can be incorporated into 
the SED for several intrinsic and extrinsic reasons, including 
(in approximate order of overall severity): 

• poor-quality raw data, e.g. saturated images or un- 
flagged cosmic rays; 

• blending, particularly among catalogues which have in- 
tegrated fluxes over different areas (e.g. the point-spread 
function of Hipparcos is a very different size to that of 
IRAS), affecting binary and multiple stars, and objects 
which are in the same line of sight as background objects 
with very red colours (e.g. redshifted galaxies in infrared 
surveys); 

• poor background subtraction, which mainly affects in- 
frared observations of sources at low Galactic latitudes or 
in other regions of nebulosity; 

• intrinsic variability among stars; 

• inaccurate source matching due to large proper mo- 
tions; 

• the accuracy of photometric corrections to each survey's 
base photometric system. 

The issues become particularly problematic when compar- 
ing optical and mid-infrared catalogues: in the mid-infrared 
stars typically present fainter detections compared to higher 
backgrounds, and beam sizes are typically larger (hence in- 
clude more objects) . Additional errors come from the model 
parameters which arise from our assumption that stars have 
solar metallicity, a given stellar mass and no interstellar ex- 
tinction. Bearing this in mind, we have assumed that each 
photometric point has an arbitrary absolute error of 10% 
and computed a x 2 -based goodness-of-fit measure based on 
this uncertainty. We also do not give error estimates for our 
temperatures and luminosities, as they would essentially be 
meaningless. 

The bad data in our combined multi-wavelength cata- 
logue is mostly of sources which have saturated in the var- 
ious input survey catalogues. The WISE catalogue, in par- 
ticular, suffers from saturation. Some photometry flagged 
as good does not match detections in other bands (e.g. 
HIP 24 436 has W 2 = 1.987 ± 0.009 mag, while Wi,s,4 = 
0.739, 0.001 and -0.030 mag). Conversely, some photometry 
flagged as uncertain does not decrease the goodness-of-fit 
of the model SEDs and is therefore sufficiently accurate for 
the purpose of identifying infrared excess. We have so far 
included all WISE data, regardless of its uncertainty and 
we must now remove the points we believe to be in error. 

Unfortunately, this is a particularly recalcitrant dataset 
to remove bad data from: we wish to keep points which fit 
badly due to intrinsic variability (as these will, on average, 
cancel out across the SED), but remove datapoints which 
have incorrect fluxes. We have opted to apply a number of 
sequential cuts to remove bad quality data. Since altering 



one band affects the model fit of the others, we must care- 
fully design these cuts to minimise errant removal of good 
quality data. At each step, we have visually examined the 
SEDs of the objects with the worst-fitting models and de- 
vised a cut which removes the dominant contribution of bad 
data. The removed data was examined and the cut applied 
if it did not remove any plausibly accurate data. The details 
of these cuts are given in the Appendix. 

While these cuts have not removed every single bad data 
point, they have removed the vast majority of bad data, 
providing a much cleaner data set with which we can work. 
This has sadly meant removing stars where there was not 
sufficient data to provide a robust fit, meaning the origi- 
nal Hipparcos catalogue has been reduced from 117 956 to 
107 619 objects. 



2.4 Denning infrared excess 

Now that we have removed bad data from our catalogue, 
we can calculate the amount of infrared excess present for 
each star. Having already performed SED fitting, creating 
a measure of infrared excess becomes a trivial exercise in 
comparing the fitted model with observations. Providing the 
best metric(s) to quantify infrared excess is more difficult. 
We adopt two techniques. 

Our first metric simply takes the ratio of observed 
to SED-modelled flux of all the data longward of 2.2 /im 
(WISE, MSX, AKARI and IRAS), and averages them 
together. This provides a single number (-Eir) that describes 
the average excess in the 3-25 /jm region, relative to the un- 
derlying photospheric model, which can be described math- 
ematically as: 

rpobs / rpmodcl 

Em= £ , (2) 

A>2.2fim 

where F° bs and F™ odel are the observed and modelled fluxes 
at frequency v (corresponding to wavelength A) and n b s is 
the number of observations at wavelengths >2.2 /im. 

Our second metric assumes that the infrared excess is 
due to reprocessed stellar light (i.e. ignoring background 
infrared emission and foreground circumstellar extinction). 
We approximate the amount of reprocessed light as being 
the integrated observed flujfl longward of 2.2 /xm, minus the 
integrated model flux over the same region. We can take 
this as a fraction of the underlying stellar flux, under the 
assumption that the total energy output (in Jy) of the star 
is not affected by circumstellar reprocessing of light. Math- 
ematically, we can then define the fraction of stellar light 
reprocessed into the infrared (Lir/L*) as: 

_ JZ»J F ° bs ~ F r dcl ) ^ 

We also define a wavelength, A poa k, as where vF° hs — vF™° Ael 
reaches a peak. The precision with which we can define A pca k 
depends strongly on the amount of available data. 



2 Here, the flux is denned in energy terms, i.e. J F v Au or 
J \F\d\. 
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Figure 1. Density-coded Hertzsprung-Russell diagram for the clipped Hipparcos dataset, based on the BT-Settl models. Top panel: 
stars with Hipparcos distances of <1 kpc. Bottom panel: stars with Hipparcos distances of <300 pc with parallax errors of <30% and 
photometric data at >2.2 (im. Darker/redder colours show regions with a greater number of stars. 
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3 MASTER CATALOGUE AND H-R 
DIAGRAM 

3.1 Presenting the catalogue 

We are now in a position to list our data in a master cata- 
logue containing the fundamental parameters of the Hippar- 
cos stars. We do so in Table [2] which lists: 

• Column 1: the Hipparcos identifier for the object; 

• Columns 2 & 3: the co-ordainates of the object; 

• Columns 4 & 5: the parallax distance to the object and 
its associated fractional error; 

• Column 6: the modelled effective temperature of the 
object; 

• Column 7: the modelled luminosity of the object; 

• Columns 8 & 9: the shortest and longest wavelength for 
which we have data; 

• Columns 10-30: the ratio of the observed to modelled 
flux at each wavelength, such that unity represents a per- 
fect match to the model; 

• Column 31: nm,, the number of observations at A > 
2.2/xm; 

• Column 32: -Eir, the average excess at those wave- 
lengths, as defined above; 

• Column 33: Lir/L*, the fraction of the object's lumi- 
nosity reprocessed into the infrared, as defined above; 

• Column 34: A pca k, the wavelength at which the infrared 
excess reaches a peak, as defined above. 



The catalogue is displayed as a Hertzsprung-Russell 
(H-R) diagram in Figure [T] The top panel of this Figure 
shows the Hipparcos sample to separate out quite cleanly 
into the two traditional populations: the main sequence 
stars, which are largely complete for stars brighter than a 
few solar luminosities, and the giant branch stars, which lie 
to cooler temperatures. The concentration of stars on the 
giant branch is due to two groups of stars. The first be- 
ing the horizontal branch stars (which, since the sample is 
largely metal-rich, form a red clump). The second being the 
RG B bump, the posi tion which is also affected by metallic- 
ity (|Cho fc Leell2002l ). A significant scatter of stars is seen 
away from these two groups, which is not necessarily real 
and which we discuss in Section T3. 2 1 

3.2 Distance-limiting the sample and associated 
biases 

Despite the removal of a significant amount of bad data, 
there is still a large amount of scatter in the H-R diagram. 
There are four reasons for this. Firstly, a large number of 
Hipparcos stars have relatively poor-quality datc0. These 
are mostly stars where there is insufficient infrared data con- 
straining the SED (these can be identified as those stars in 
Table [5] where niR, the number of measurements at A > 2.2 
fim, is small). This leads to a large fraction of the scatter 
observed in the H-R diagram, and to vertical artefacts (con- 
centrations and rarefactions) on the BT-Settl model grid 
spacing. 

3 We refer here to noise over the entire SED, rather than one or 
two clearly mismatching 'bad' points. 



Secondly, stars which do have infrared excess or suf- 
fer from substantial interstellar extinction scatter towards 
cooler temperatures (and, in the case of interstellar extinc- 
tion, lower luminosities) in the H-R diagram, as their optical 
light is either reprocessed into the infrared or scattered out 
of the line of sight. 

Thirdly, and perhaps most importantly, scatter arises 
from uncertainty in the Hipparcos parallax, which smears 
objects vertically in the H-R diagram. T his also leads to 
the Lutz-Kelker bias l|Lutz fc Kclkcr 197jl). This bias occurs 
due to the inversion of parallactic angle to obtain a distance. 
As the measurement error is in parallax, this preferentially 
scatters objects to smaller distances. Hipparcos data suffers 
from this significantly. We use the benchmark of 17.5% error 
in parallax (which correspond s to an average 30% de viation 
in a set of stellar luminosities; lOudmaiier et alii 19981 ) as our 
figure of merit. Of the 107 619 stars in our final sample, only 
49 188 have parallax errors less than this value. This rather 
severe limitation reduces the usefulness of the sample in ex- 
amining stellar populations, particularly for the relatively- 
rare stars on the upper giant branches. For many applica- 
tions a wider sample, with increased Lutz-Kelker bias is 
preferable. We therefore continue to include objects suscept- 
able to significant Lutz-Kelker bias, but warn the reader to 
be mindful of its existence. 

Finally, an additional distance error is present in red 
giants, where changes in brightness across the stellar sur- 
face (which covers a finite solid angle) cause a measurable 
astrometric shift. This can be misinterpreted as a parallac- 
tic shift, leading to much smaller Hipparcos distances than 
their true distance. This is perhaps best observed in the case 
of W Lyn, which has a Hipparcos par allax of 21.53 ± 8.06 
mas, despite being several kpc distant l|lta et al.ll200ll ). This 
is an extreme case, though we warn the reader that no sup- 
posedly volume-limited sample of any consequence will be 
clean of all intruding objects for these reasons. 

For the remainder of our analysis we adopt two volume- 
limited samples, which are subject to these biases at differing 
levels. The first is limited to stars with parallax distances of 
<300 pc with <30% parallax errors, which have a wide range 
of data which cover the SED well (A max > 2200«m). This 
subset of data still contains 46 869 of the original Hipparcos 
stars, of which 34 660 have parallax errors below 17.5%. The 
second sample is distance limited to <200 pc with <30% 
parallax errors, with the same requirement that in which 
32 741 out of 33 898 stars have errors <17.5% in parallax. 
We refer to these in the discussion as the 300-pc and 200-pc 
samples. 



4 DISCUSSION 

4.1 Stellar isochrones 

Much can be made of the H-R diagram in terms of the star 
formation history of the local neighbourhood. However, to 
do so thoroughly requires a more in-depth analysis than we 
are able to provide in this work. As a cursor y analysis, we 
present solar-metallicity Padova isochrones (jMarigo et al.l 
120081 ; iBertelli et ail 120081 ) for a variety of ages in Figure [2] 
The H-R diagram is well-described by a population of mixed 
age, but of near-solar metallicity. 
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Table 2. Fundamental parameters and infrared excess for Hipparcos stars. The columns are described in the text. The complete table 
is available in the online version of the paper. 



HIP RA Dec 


d 


8d/d 


T e s 


L 


Coverage 


(nm) 


SDSS u ■ 


■IRAS [25] 


nm 


Sir. Liu 


-^pcak 


(J2000) (J2000) 


(PC) 




(K) 


(L©) 


Start 


End 


Excess 


■ Excess 




L, 


(H 


1 0.00091 +01.08901 


219.78 


0.29 


6400 


8.73 


354 


2200 


0.872 


■ ■ ■ 0.000 





0.000 





2 0.00380 -19.49884 


47.96 


0.05 


3300 


3.11 


354 


623 


1.077 


■ • ■ 0.000 





0.000 





3 0.00501 +38.85929 


442.48 


0.15 


8968 


374.87 


420 


8610 


0.997 


■ • ■ 0.000 


1 


2.300 0.0019 


8.6 




20000 10000 5000 3000 

Temperature (K) 

Figure 2. Density-coded Hcrt zsprung-Russell diagram for the 200-pc sample (greyscale). Overplotted are solar-metallicity isochrones 
from the Padova models llMarigo et alj|2008t iBertelli et al.ll2008l) at 10, 20, 30 and 50 Myr (solid, red lines); 100, 200, 300 and 500 Myr 
(long-dashed green lines); 1, 2, 3 and 5 Gyr (short-dashed blue lines); and 10 Gyr (dotted magenta line). The thin red line to the left of 
the main sequence is a zero-age isochrone at [Fe/H] = -1 to illustrate the blueward shift caused by decreasing metallicity. Black arrows 
show the effect of dereddening individual sources by E(B — V) = 0.1 mag. 



Interstellar reddening does not appear to be a great 
cause of concern in the 200-pc sample in general, with the 
majority of stars lying within the bounds of the isochrones 
with only a few tenths of a magnitude of dereddening at 
most. We show on Figure [2] the effect that a reddening of 
E{B — V) — 0.1 mag has on particular sources chosen at 
a variety of different temperatures. We can see here that 
cooler sources are largely unaffected by this modest redden- 
ing, but that the effect becomes much more severe as the 
peak of the SED becomes bluer and the short-wavelength 
photometry available fails to constrain the SED. This may 
lead to underestimates of the temperatures and luminosities 
of some of the hotter stars. Particularly, errors can be large 
if short- wavelength data is unavailable or unusable. 



4.2 Comparison to spectroscopic temperatures 

With this in mind, we can check the consistency of our re- 
sults, by comparing them to the s pectroscopically-derive d 
temperatures of the NStars project (IGrav et al.ll2003l . 1200^ 1. 
This project identifies the basic parameters of stars within 
40 pc of the Sun by fitting moderately-high-resolution spec- 
tra, thus the results are unbiased by interstellar reddening. 
Figure [3] shows the ratio of their temperatures to ours for 
407 stars we have in common. Examination of the three 
outliers (HIP 35 550, 59 199 and 71 957) in this figure show 
obvious problems with the 2MASS photometry that were 
missed by our bad data cuts. Neglecting these, the average 
temperature is consistent to <0.22% (i.e. the error on the 
mean). The standard deviation of results is 4.4% overall and 
decreases slightly toward lower temperatures. 

Interstellar reddening should not affect the spectro- 
scopic temperature, but imparts an apparent cooling to 
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NStars spectroscopic temperature (K) 

Figure 3. Comparison of our temperatures with those derived 
spectroscopically from the NStars project. Lines show the ap- 
proximate deviations expected for E(B — V) ~ 0.0, 0.1 and 0.2 
mag. 



the photometric temperature, scattering points above unity. 
Four stars have a spectroscopic temperature >15% higher 
than the photometric temperature: HIP 71 193, 84 379 (8 
Her), 93 805 (A Aql) and 98 495 (e Pav). The latter two scat- 
ter in this direction due to excess flux in the WISE 4.6-fim 
band. The former two have SEDs that are poorly constrained 
by the input photometry. Unsurprisingly, we therefore find 
negligible reddening among stars within 40 pc. 

In principle, with sufficiently-good-quality photometry, 
one could compute a three-dimensional extinction map of 
the Solar Neighbourhood by comparing spectroscopically- 
derived temperatures to photometrically-derived tempera- 
tures. The lack of self-consistent photometry for the Hip- 
parcos stars probably prevents such a determination here, 
but may become possible in the Gaia era. 



4.3 Grouping objects by type 

The AKARI-Hipparcos catalogue of llta et al.l (|201fJ ) mea- 
sures excess at 9 and 18 /jm, and contains grouping infor- 
mation for 2787 stars commonly exhibiting infrared excess, 
including carbon stars, red giants, supergiants, S-type stars, 
etc. AKARI only detected 44% of the Hipparcos stars, so 
can be regarded as a selected subset of the Hipparcos sam- 
ple, subject to its own biases. We remind the reader that the 
Hipparcos sample does not include all of the local (<300 pc) 
optically-obscured giant stars, and that the compilation of 
stellar types listed in Ita et al. are only as complete as the 
literature from which they are based. The catalogue from 
llta et all l|2010h cannot therefore be treated as a definitive, 
complete list of each type of source, nor does it purport 
to be such. Of the 2787 classified stars in the list of Ita et 
al., we retain 2764 after removing bad data. Of these, only 
749 meet the criteria for our 300-pc sample and 293 for our 

200-pc sample. 

The objects classified in llta et al.l (|2010t ) within the 300- 
pc sample are presented on a H-R diagram in Figure [4] In 
general, the different types of stars match up well with their 
expected locations. M giants and S stars lie on the upper 
giant branch, along with the majority of carbon stars (we 
have not made an effort here to separate intrinsic from ex- 
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Figure 5. Top panel: Hertzsprung-Russell diagram of the 300-pc 
Hippacros sample. Bins are colour-coded by infrared excess (Em), 
with the colour-coding denoting the object with the greatest ex- 
cess in that bin: light blue indicates no excess (-Eir, = 1), green 
indicates moderate excess (-Em ~ 2), red strong excess (Ejr ~ 5) 
and black extreme excess (capped at Em = 10). Middle panel: as 
top panel, but showing the Galactic distribution of those sources. 
Bottom panel: as top panel, removing sources with |b| < 5° and 
within 5° of the IRAS missing strip. 



trinsic carbon stars (|Van Eck et al.|["l998h due to the incom- 
plete nature of any determination). Be stars (which include 
a variety of hot, emission- line stars) are located in the upper- 
left of the H-R diagram, but scatter toward temperatures 
cooler than the main sequence due to reprocessing of stel- 
lar light into the infrare d by the circumstellar excretion disc 
(cf. iKastner et al]|2006h . Pre-main-sequence (pre-MS) stars 
likewise mostly lie toward the cooler side of the main se- 
quence due to reprocessing of their optical emission into the 
infrared, lowering the effective temperature. 

The top panel of Figure[5]shows a similar H-R diagram, 
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Figure 4. Density -coded He r tzspr ung-Russell diagram for the 300-pc sample (greyscale, from bottom panel of Figure[TJ. The ovcrplotted 
symbols are from llta et aL and show known examples of the following. Black filled circles: pre-main-sequence stars — cyan 

downward-pointing triangles: Be stars — large, magenta upward-pointing triangle: Wolf— Rayet star — blue crosses: M-type (super- 
giants — green asterisks: S-type giants — red filled squares: carbon stars. 



indicating the location of the stars with greatest infrared ex- 
cess, as measured in terms of average fractional excess over 
the mid-infrared spectrum (Em). The two regions of high 
stellar density (the main sequence above 6000 K, and the 
red clump) contain the largest number of infrared-excessive 
stars. These stars are not usually truly excessive, except the 
case of several pre-MS stars, but are instead artifically red- 
dened or suffer from source confusion in the infrared. The 
majority of these stars are located near the Galactic Plane, 
or lie in or near the IRAS missing stripe and thus suffer from 
poor infrared photometry. Some stars in the Plane may suf- 
fer from sufficient interstellar extinction to appear to have 
infrared excess, even at distances as small as 300 pc. Remov- 
ing sources within 5° of the Galactic Plane and the IRAS 
missing strip yields the H-R diagram at the bottom of Fig- 
ure [5] The number of highly-excessive sources is greatly re- 
duced, with the few remaining sources of high excess located 
predominantly in the Gould Belt. 

More pertinently, two further groups of infrared- 
excessive stars in Figure \S\ are largely unchanged by this 
process. The first, above the main sequence a t arou nd 10 000 
K, are identified as Be/Ae stars bv llta et alj l|2010l ): the sec- 
ond are the cool, luminous stars near to and above the RGB 
tip. This second group are ident ified as M giants and S and 
carbon stars by llta et al.l l|2010l ). and are likely to entirely 
be dust-producing AGB stars. 

In Fig ure [6] we show the different types of identified 
stars from llta et all (|2010h as a function of our two mea- 



sures of infrared excess: Em and Lm/L*- While Lm/L* is a 
more physical measure of infrared excess, Em is clearly more 
effective at separating out the infrared excessive stars from 
the bulk of the population, particularly for the Be stars. 

Figure [S] also shows a signficiant scatter below _Eir, = 1. 
This represents a supposed deficit of infrared flux in these 
stars, and tends to be much more prevelant in the cooler 
stars. This is largely due to decreased sensitivity in stars fur- 
ther down the main sequence, but some upper-giant-branch 
stars also have infrared deficits. In this case, scatter can 
be introduced by stellar variability, as photometry is not 
usually averaged or taken contemporaneously, and does not 
imply either an instantaneous or time-averaged deficit in in- 
frared flux. 

4.4 Giant stars with excess 



Figure [7] shows that scatter on the giant branch gen- 
erally decreases as we approach the AGB tip, due to in- 
creased sensitivity on brighter sources. We can also see the 
substantial number of sources which have scattered to lower 
temperatures and higher values of Em (also Figure [6}, indi- 
cating reprocessing of optical light be circumstellar dust. 
The amount of excess around these stars can be corre- 
lated with their dust-production rate and hence their mass- 
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Table 3. Literature spectroscopic and variability information for luminous (>850 Lq) giant stars with detected circumstellar emission. 
Details of columns and explanations of acronymns are listed in the text. 



HIP Name Spectral Temper- Lumin- Variability LRS SWS E m 

Type ature osity Type Amplitude Band Source Period Source Classif- Classif- 









(K) 


(Lq) 




(mag) 






(days) 




ication 


ication 




1728 


T Cet 


M5-6SIIe 


3329 


7442 


SRV 


1.9 


V 


G 


159 


G 






2.45 


8565 


TT Per 


M5II-III 


3228 


1579 


SRV 


1.1 


P 


G 


82 


G 


SE7 




1.89 


9234 


V370 And 


M7III 


2948 


3831 


SRV 


1.01 


H 


G 


228 


G 






4.44 


13064 


Z Eri 


M5III 


3354 


2334 


SRV 


1.63 


V 


G 


80 


G 


SE8 




1.79 


17881 


SS Cep 


M5III 


3158 


5104 


SRV 


1.1 


P 


G 


90 


G 


SE3 




1.53 


24169 


RX Lep 


M7 


3256 


3764 


SRV 


2.4 


V 


G 


60 


G 


SE6 




2.19 


25194 


SW Col 


M1III 


3661 


921 


V 


0.34 


V 


G 






SE8 




2.40 


27989 


a Ori 


M2Iab: 


3659 


73524 


SRV 


1.3 


V 


G 


2335 


G 




2.SEcp 


2.88 


28166 


BQ Ori 


M5IIIv 


3192 


917 


SRV 


2.1 


V 


G 


110 


G 






1.95 


28816 


SS Lep 


A1V+M6II 


4347 


2672 


Bin 


0.24 


V 


G 










12.60 


28874 


S Lep 


M5III 


3187 


3415 


SRV 


1.58 


V 


G 


89 


G 


SE6t 




3.79 


36288 


YLyn 


M6SIb-II 


3110 


5249 


SRV 


2.5 


V 


G 


110 


G 






2.10 


38834 


V341 Car 


M0III 


3326 


1580 


V 


0.9 


V 


G 





G 






5.88 


41037 


e Car 


K3:IIIv+? 


4209 


14086 


Bin 


0.08 


H 


G 





G 






1.34 


42489 


RV Hya 


M5II 


3200 


1884 


SRV 


1.3 


V 


G 


116 


G 






2.36 


43215 


AK Pyx 


M5III 


3410 


1499 


V 


0.42 


H 


G 





G 






1.77 


44050 


RT Cnc 


M5III 


3192 


2225 


SRV 


1.48 


V 


G 


60 


G 


SE3 




2.17 


44862 


CW Cnc 


M6 


2909 


2228 


V 


1.2 


P 


G 






SE3 




2.93 


45058 


RS Cnc 


M6IIIase 


3122 


5282 


SRV 


1.5 


P 


G 


120 


G 






2.44 


46806 


R Car 


M6.5IIIpev 


2800 


4164 


Mira 


6.6 


V 


G 


309 


G 


SE1 




2.38 


48036 


R Leo 


M8IIIe 


1995 


1493 


Mira 


6.9 


V 


G 


310 


G 


SE2 




2.43 


51821 


U Ant 


C5,3(Nb) 


3317 


5819 


V 


0.9 


p 


G 






SiC+: 




2.29 


52009 


U Hya 


C6.5,3(N2) 


3400 


3893 


SRV 


2.4 


B 


G 


450 


G 


SiC 




2.22 


53809 


R Crt 


M7III 


2491 


8591 


SRV 


1.4 


P 


G 


160 


G 


SE3t 




3.79 


57607 


V919 Cen 


M7III 


3094 


7766 


SRV 


0.58 


H 


G 










3.17 


61022 


BK Vir 


M7III: 


2889 


2706 


SRV 


1.52 


V 


G 


150 


G 


SE4t 




2.64 


63642 


RT Vir 


M8III 


2602 


1804 


SRV 


1.29 


V 


G 


155 


G 


SE3t 


2.SEa 


4.12 


64569 


SW Vir 


M7III 


2918 


4917 


SRV 


1.5 


V 


G 


150 


G 


SE3t 




3.22 


68357 


RW CVn 


M7III: 


3141 


973 


SRV 


1.1 


P 


G 


100 


G 


SE2: 




2.98 


68815 


9 Aps 


M6.5III: 


3151 


3879 


SRV 


2.2 


P 


G 


119 


G 


SE5t 


2.SEb 


3.48 


69816 


U UMi 


M6e 


3018 


1821 


Mira 


5.9 


V 


G 


331 


G 


SE2 




2.20 


70401 


RX Boo 


M7.5 


2581 


8196 


SRV 


2.67 


V 


G 


162 


G 


SE3t 


2.SEa 


3.02 


70969 


Y Cen 


M7III 


2907 


5317 


Irr 


1.1 


p 


G 


180 


G 


SElt 




1.61 


71802 


RW Boo 


M5III: 


3148 


3010 


SRV 


1.5 


V 


G 


209 


G 


SE7 


2.SEb 


1.77 


72208 


EK Boo 


M5III 


3333 


5587 


SRV 


0.38 


H 


G 










1.46 


76423 


T4 Ser 


M5II-III 


3165 


5264 


SRV 


1.18 


V 


G 


100 


G 


SE4 




1.88 


77619 


ST Her 


M6-7IIIaS 


3071 


6270 


SRV 


1.5 


V 


G 


118 


G 


SE1 


2.SEa 


2.70 


78574 


X Her 


M8 


3152 


2765 


SRV 


1.1 


V 


G 


95 


G 


SE6t 


2.SEb 


5.65 


80488 


U Her 


M7III 


2700 


4438 


Mira 


7 


V 


G 


406 


G 


SE4 


2.SEc 


4.16 


80704 


g Her 


M6III 


3261 


4056 


SRV 


2 


V 


G 


89 


G 






1.46 


84345 


a Her 


M5Iab: 


3351 


15368 


Bin 


1.26 


V 


G 








l.NOp 


1.45 


86527 


BM Sco 


K2.5Iab: 


3676 


949 


SRV 


1.9 


P 


G 


815 


G 






5.69 


94162 


SZ Dra 


M 


3173 


923 


V 


1 


P 


G 






SE5 




3.21 


95413 


CH Cyg 


M7IIIv 


2687 


4316 


Bin 


2.89 


V 


G 








2.SEc 


5.39 


95902 


AF Cyg 


M4 


3305 


1646 


SRV 


2 


V 


G 


93 


G 


SE3t 




2.29 


98031 


S Pav 


M7IIe 


2752 


5563 


SRV 


3.8 


V 


G 


381 


G 


SE2t 


2.SEa 


2.28 


99082 


V1943 Sgr 


M7III 


2752 


5813 


V 


2 


p 


G 






SE2t 


2.SEa 


2.04 


99990 


RT Cap 


C6,4(N3) 


3245 


2683 


SRV 


2.8 


p 


G 


393 


G 


SiC+: 




1.98 


100935 


T Mic 


M7III 


2856 


7708 


SRV 


1.9 


p 


G 


347 


G 


SElt 


2.SEa 


2.46 


101810 


EU Del 


M6III 


3227 


1585 


SRV 


1.11 


V 


G 


60 


G 


N 




1.77 


104451 


T Cep 


M5-9 


2866 


6767 


Mira 


6.1 


V 


G 


388 


G 


SE1 


2.SEa 


2.85 


107516 


EP Aqr 


M8IIIv 


3056 


2651 


SRV 


0.45 


V 


G 


55 


G 


SE5t 


2.SEb 


4.27 


108928 


TW Peg 


M7.5IIIv 


3145 


5027 


SRV 


0.9 


p 


G 


929 


G 


SE6t 




6.11 


110396 


DZ Aqr 


M 


3055 


2454 


V 


1.1 


V 


G 








2.SEb 


3.44 


110428 


BW Oct 


M7III 


2849 


3592 


V 


0.9 


p 


G 






SE5t 




2.32 


114318 


Y Scl 


M6III 


3039 


872 


SRV 


1.6 


p 


G 






SE7 




2.48 


114404 


V345 Peg 


M3 


3345 


2795 


V 


0.37 


H 


G 










1.72 


117245 


TX Psc 


C7,2(N0) 


3451 


5693 


V 


0.41 


V 


G 






N 


l.NC 


2.28 


118188 


R Cas 


M7IIIc 


2187 


2219 


Mira 


8.8 


V 


G 


431 


G 


SE5t 


2.SEb 


4.73 
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Table 4. Literature spectroscopic and variability information for less-luminous (<1000 Lq) giant stars with detected circumstellar 
emission. Details of columns and explanations of acronymns are listed in the text. 



HIP Name Spectral Temper- Lumin- Variability LRS SWS E m 

Type ature osity Type Amplitude Band Source Period Source Classif- Classif- 









(K) 


(L Q ) 


(mag) 






(days) 




ication ication 


893 


AC Cet 


M3III 


3413 


793 


SRV 0.33 


V 


G 






2.71 


39751 


RU Pup 


C5,4(N3) 


3323 


394 


SRV 1.9 


P 


G 


425 


G 


2.92 


43438 


RS Cam 


M4III 


3298 


739 


SRV 1.8 


V 


G 


89 


G 


1.94 


50916 


HR 4091 


K4III 


4057 


459 


V 0.03 


H 


H 






2.83 


52656 


TZ Car 


C (R5) 


3346 


326 


SRV 1.7 


P 


G 


69 


G 


2.90 


57800 


RU Crt 


M3 


3054 


681 


V 1 


P 


G 




G 


SE3: ■ ■ ■ 2.04 


59389 


HD 105822 


K0/K1III 


4464 


234 


None 0.03 


H 


H 






1.91 


59458 


68 UMa 


K5III 


4478 


256 


None 0.04 


H 


H 






2.70 


71568 


HR 5464 


K4III 


4214 


727 


None 0.05 


H 


H 






1.68 


112155 


BD Peg 


M8 


3147 


736 


SRV 0.9 


P 


G 


78 


G 


2.56 
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Figure 6. Measures of excess for different kinds of identified stars 
(see text for definitions of i^jR and -Ltr). Symbols are as in Figure 

HI 



loss rate. By identifying and characterising individual stars 
which lie above the general scatter in Figure and includ- 
ing optically-obscured sources missed by Hipparcos, we can 
make a theoretically complete census of dust producing stars 
within 300 pc. 

While that is beyond the scope of this paper, we do re- 
port on giant stars which are observed to have significant 
excess. We define this by Em > 2.65 — log(L)/3 (see also 
Figure [7|, which is chosen to identify excesses of >2<r at all 
luminosities. We list these stars in Tables [3] (luminous stars, 
above 850 Lq) and|4] (stars below 850 Lq). In these tables, 
we have also listed common names and spectral classification 
from SIMBAD, and variability info rmation, sourced from ei- 
ther the GCVS (|Samus et al.ll2006j ; denoted 'G') or Hippar- 
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Figure 7. Excess among giant stars. Colour scale denotes mod- 
elled stellar effective temperature: red points are coolest, blue 
points are warmest. Red plus signs show stars from I GroenewegerJ 
(|20L2j) that were identified to have optical depths of ry > 10 — 5 , 
blue crosses show those stars that were not. The dotted black line 
shows our definition of those stars with infrared excess. 
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Figure 8. SED of EU Del. Large red points show literature pho- 
tometry, small blue points show literature spectra. The black 
line indicates the Rayleigh- Jeans tail expected for a dustless or 
'naked' star. References are listed in the text. 



cos (denoted 'H') catalogues. Variability periods and ampli- 
tudes (either in the Hipparcos band, denoted 'H'; Johnson 
V-band, denoted 'V; or as a photographic measurement, 
denoted 'p';) are shown where available. Variability types 
are listed as follows: Mira = Mira variable; SRV = semi- 
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regular variable; Irr = irregular variable; V = unclassified 
variable; Bin = (eclipsing) binary; None = no appreciable 
variability. The spectral classifications of any LRS and SWS 
spectra of these sources are also given ijSloan fc Price|[l99il ; 



ISloan et al. 19981. 120031; based on the classification method 



of lKraemer et al.ll2002T ). 

Table [3] contains some well-known targets, which are 
known to have substantial mass loss. The common SE clas- 
sification of the LWS and SWS spectra shows that most of 
these stars are known to have silicate features in their spec- 
tra and are therefore dust producers. Others, such as SS 
Lep, e Car and a Her, are known to be binaries, therefore 
their SEDs may not be well-repesented by a single blackbody 
(most giant stars, however, should outshine any companions 
at all wavelengths: note that a Her in particular is known to 
be dust producing; iTatebe et al.ll2007l ). Several others with- 
out LRS or SWS spectra are likely to be mass-losing stars, 
but without infrared spectra it is difficult to tell. 

One outstanding example is present, however: that of 
EU Del (HIP 101810). Figure |3] shows its SED. Liter- 
ature photometry for this SED comes from the Hippar- 
cos, Tycho, 2MASS, AKARI and IRAS catalogues already 
mentioned. Further photometry was sourced from Mermil- 
liod's catalogue of homogeneous means (UBV; VizieR on- 
line data catalogue H/168) , the Carlsberg Meridian Cata- 
logue (V;\E vans et al . 1999), the catalo gue of infrared obser- 
vations, including the Revised AFGL cata logue (4-20 /im; 
iPrice fc Murdockll 19831; ICezari et"al] 11993ft , DIRBE (2.2- 
100 mm lPrice et aJj|201Ch . and the AKARI FIS catalogue 
(65, 90 and 140 /xm; VizieR online dat a catalogue 11/298) . 
Literatu re spectroscopy for t he optical (Valdcs et al. 2004), 
J-band (I Wallace et all 120001). iJ-band dMever et al.l I1998T ) 
and if -band (|Wallace fc Hinkld 1 1997 ) are also sh o wn, a s 
is the IRAS LRS spectrum from ISloan fc Pried (|l99ST ). 
ISloan fc Price] (|l998l ) classify EU Del as a naked star, as 
it shows no silicate feature. However, it is found to have 
substantial excess in both the IRAS 12- and 25-/xm bands 
and the infrared spectrum clearly shows a ri s e abo ve a black- 
body toward longer wavelengths. IWu et all (1201 ll ) place the 
star at a metallicity of [Fe/H] = -1, making it reminiscent of 
the featureless excesses we have pr eviously found in metal- 
poor g lobular cluster giant stars l|McDonald et al.l l2010al . 
1201 lal ldT) . We have previously attributed this to metallic 
iron dust on chemical and mineralogical grounds, however it 
is spectrally indistinguishable from amorphous carbon dust 
and silicate dust compo s ed primarily of la rge grains (cf. 
iHofner fc Andersen! [20071 : iNorris et all 120121 '). EU Del may 
therefore be a unique nearby testbed in which to determine 
which dust species is causing these unexplained featureless 
excesses. 

Table [4] contains a few stars with suspected infrared ex- 
cess, which we examine more closely, in order to identify the 
luminosity at which dust formation (as traced by infrared 
excess) begins. We investigate these individually here. 



• AC Cet: Considerable excess exists in the AKARI 9- 
and 18- Atm bands , and the IRAS 12- and 25-/xm bands. 
IChen et all l|2004h note that there is another evolved star 
within the IRAS PSF, but the excess in AKARI sug- 
gests that AC Cet does indeed have circumstellar dust. 
iKwok et al.l (|l997l ) classify the source as having a class 



'C LRS spectrum, denoting a ll-/im SiC feature, how- 
ever this feature is tentative in this source, at best. 
• RU Pup: AKARI, WISE and IRAS data all show con- 
siderable excess at wavelengths longer than 4 /im. There is 
significant scatter in the optical photometry, which leads 
to a poor estimation of the temperature and luminosity 
for this star. This may be partly due to its carbon-richness 
and partly due to its variability. iBergeat fc Chevallierl 



__P< 

(|2005T ) place this star via two means at 2680 or 2875 K 
and 455 or 610 pc, which makes it considerably cooler and 
more luminous (1715 or 3649 Lq) th an we model. This 
is corroborated by its long period (cf. Ilta et al1l2004l) . It 
therefore probably suffers from the pulsation-induced dis- 
tance errors we describe in Section [2. 21 
• RS Cam: This star also shows excess in the AKARI, 
WISE and IRAS data. Its short p eriod suggests its lu- 
minosity is correctly determined (cf. Ilta et aU 2004 ). The 



LRS spectrum shows weak silicate emission (|Kwok et alj 
Il997t ). 

• HR 4091: This source is modelled using DIRBE and 
IRAS data only. It shows marginal excess between 4 and 
21 /im, and considerable excess in the 25-^m IRAS band. 
It is at very low Galactic latitude (b = —0.5°) and thus 
suffers from considerable contamination from surrounding 
sources. We therefore do not believe this excess is real. 

• TZ Car: The 8- to 25-/im data for this star show con- 
siderable exc e ss. R eprocessing of the Hipparcos data by 
iKnapp et al.l (|200lh suggest the distance for this star is 
roughly correct. At a Galactic latitude of b — —5.8°, TZ 
Car may suffer from some extinction, but it is likely that 
the excess and parameters are sufficiently correct to say 
that this star is losing mass. 

• RU Crt: A known mass-losing star, this star shows mod- 
erate excess in the AKARI bands and substantially more 
excess in the IRAS bands. At 132 pc, it is possible that 
more-extended emission is missed by AKARI: sources up 
to roughly this distance may have some emission outwith 
the AKARI beam at 25 /im (see Section [4. 5p . 

• HD 105 822: This star is in a region of high projected 
stellar density (b = —5.7°). The amount of excess for 
this star is inconsistent across the infrared data, vary- 
ing among the surveys and bands. A dubious J-band 
flux probably suggests more infrared excess than is truly 
present. The apparent infrared excess in this case is likely 
due to source blending and confusion, coupled with poor- 
quality short-wavelength photometry. 

• 68 UMa: This star is mistakenly classified as excessive 
due to a mismatch between the Hipparcos /Tycho magni- 
tudes and those from the SDSS, which differ by approxi- 
mately a magnitude, despite no variability being detected 
by Hipparcos. This has lead to an excess being determined 
incorrectly. By using a variety of combinations of pho- 
tometry, we estimate that a correct temperature and lu- 
minosity of around 4000 K and 400 L© would me more 
appropriate for this source, and that it has no substantial 
reddening. 

• HR 5464: The determination of infrared excess for this 
star is based solely on IRAS data, in which excess is rela- 
tively weak (86% at 12 pm, 50% at 25 fim). DIRBE data 
suggest there is little or no infrared ex cess for this s ource . 

• BD Peg: A known mass-losing star, IKwok et al.1 (Il997l ) 
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note that sil i cate e mission is present in this object and 
iGezari et all (| 19931 ) confirm its infrared excess. 

With the exception of the carbon star TZ Car, we there- 
fore find no detectable dust production by any object below 
the luminosity of RU Crt (681 Lq). We therefore conclude 
that this represents the luminosity at which dust production 
by AGB stars begins in earnest in the local neighbourhood. 
This corroborates very well with the «700 we have pre- 
viously found in Galactic globular clusters dMcDonald et aD 
2009: lBover et al.ll2009al ; iMcDonald et al.ll2011bl ldh. 



4.5 Comparison to Groenewegen (2012) 

We n ow turn our attention to the work of I Groenewegen! 
(2012). This paper identifies several low-luminosity Hippar- 
cos stars with infrared excess which Groenewegen attributes 
to weak dust emission. These stars have much lower lu- 
mionsities (50-350 Lq) than those we find excess around, 
as Groenewegen was examining RGB stars for dust excesses 
much smaller than we have deemed accurately determinable 
in this work. 

The stars from I Groenewegen! |2012l ) which have sur- 
vived our data quality cuts are shown on Figure [7] Groe- 
newegen purposely targetted stars with low (V — I) colours, 
thus his sample does not include stars from with high values 
of Em,. It is notable, however, that all 52 stars common to 
our datasets lie within the scatter of points with no unusual 
infrared excess. Also, Groenewegen's dusty stars appear to 
have no more infrared excess than his dustless stars. 

The reasons behind this are not immediately obvious, 
but can be understood by examining the subtle differences 
between our analyses. The most striking of these is the 
choice of input data. We have included data from WISE, 
AKARI (9 and 18 /im) and IRAS (12 and 25 urn). Groe- 
newegen includes data from AKARI (including the far- 
infrared 60- and 90-/im FIS bands) and IRAS (including 
the far-infrared 60- and 100-/im bands). 

The choice of whether or not to include the far-infrared 
data is a balance of gaining sensitivity to cold dust and ac- 
quiring systematic errors due to contamination in the line 
of sight. The issue with the AKARI FIS and IRAS data 
we have excluded is the beam size (37" and 39" for AKARI 
FIS 60 and 100 /im, 

respective!^; 4 .5' x 0.7' for IRAS 12/zm; 
iMiville-Deschenes k, Lagachel 2005). These are much larger 
than the beam sizes for the other infrared data (5.6" for 
AKARI IRC 8.6 /iirQ 7.36" x 6.08" for WISE 11.6/im0) 
and much larger than the optical photometry (typically 1— 
2"). A large point-spread function full- width half-maximum 
(PSF FWHM) means there is a substantial issue with con- 
tamination from unrelated sources in the line of sight, from 
diffuse background emission, or from ISM headed by the 
star. Equally, if a giant star has a spatially-extended wind, 
this may be missed by only considering data with a small 
PSF. Groenewegen has been careful to exclude sources with 



4 From the AKARI FIS data user manual, version 1.3: 
http: //www. ir . isas .jaxa.jp/ASTRQ-F/Dbservation/IDUM/FIS_IDUM 
From the AKARI NIR data user manual, version 1.3: 
http : //www . sciops . esa . int/SA/ASTROF/docs/IRC_IDUM_l . 3 . pdf 
u From the WISE preliminary release explanatory supplement: 
http : / /wise2 . ipac . caltech . edu/ docs/release/prelim/expsup/sec4. 



strong cirrus contamination, which is the main contaminant 
in the IRAS photometry. However, his exclusions are based 
on the 100- /im images, whereas ecliptic dust is a greater con- 
taminant at 25 and 60 /im. These sources would therefore 

not be identified as contaminated. 

On examining the individual stars which IGroenewegenl 

|2012l ) claims are dusty, we find that the comparatively-large 
PSF size of IRAS and (in some cases) AKARI FIS appears 
the primary cause of the difference between our datasets. 
Typically, the IRAS 12- and 25-/im flux is systematically 
in excess of the modelled stellar photosphere compared to 
the smaller-PSF AKARI IRC and WISE photometry. As an 
example, we model HIP 44 126 (FZ Cnc) to have moderate 
excess in IRAS (52% and 35% at 12 and 25 /im), but little 
excess in WISE (8% and 6% at 11.6 and 22.1 /im). By only 
taking the IRAS data, Groenewegen naturally models this 
star as having reasonable infrared excess. In this particular 
instance, the contaminating source can be clearly identified 
as poorly-subtracted emission from warm dust in the ecliptic 
plane in the original IRAS photometry. The same is true of 
HIP 53 449, though here is the AKARI FIS data that suffers 
from contamination from the ecliptic. 

Not all of Groenewegen's dusty sources can be explained 
so easily, however. HIP 67 605 and 67 665 (AW CVn) are 
both identified as dusty by IGroenewegenl ((2012) . They lie 
quite close to each other (15' apart) but are resolved in the 
IRAS images by several beam widths. They do not suf- 
fer from substantial contamination. They are covered by 
AKARI and IRAS, but not the WISE preliminary cata- 
logue. Both sources have excess at IRAS 12 and 25 /im but 
not AKARI 9 and 18 /im. Like the majority of Groenewe- 
gen's targets, these stars lie at around 200 pc. At this dis- 
tance the AKARI 8.6-/im PSF has a FWHM of 1120 AU, 
or 2000-4000 stellar radii. Assuming the dust temperature 
approximately follows a T 4 oc R 2 law, and a stellar temper- 
ature of ~3700 K, this implies that dust falling within the 
IRAS beam but outwith the AKARI IRC beam should emit 
with a peak wavelength of A 3> 35 /im. It should therefore 
not emit significantly at 8-18 /im to cause the discrepency 
between the AKARI NIR and IRAS 12-/im fluxes. 

Nevertheless, the infrared excess Groenewegen finds 
may still be real, and still be related to a wind eminat- 
ing from the star. Two situations may cause this. The first 
case is that a cooler, detached shell surrounds the star 
and e mits only at longe r (60-100 /im) wavelengths (cf. Y 
CVn; iLibert et al. 2007). In the second case, the emission 
would be produced not by the star, but instead interstel- 
lar dust swept up in a bow shock around the astropause 
|Wareing et all 20071'). This has been seen in other n earby 
AGB stars (|Ueta et all I2006I . l2010l : lLadial et all |2010| ) and 
could be the source of the excess emission at longer (60-100 
/im) wavelengths that Groenewegen finds in several cases. 

We also acknowledge that variability may also play a 
role in this analysis. Neither Groenewegen's nor our deter- 
minations o f infrared exce s s take into account the variability 
of stars. As IGroenewegenl (|2012l ) uses some different optical 
data to us, we may find that some stars in both studies 
^a^ppear to have infrared excess simply because their optical 
piiotornetry was observed when the star was at photometric 
minimum. Conversely, excess may be missed if observations 
wer e carried out a t phot ometric maximum. 
.5c . h filroenewegen (2012) and this work probe subtly dif- 
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ferent datasets with subtly different techniques. It should 
therefore not be surprising that we find different results, 
though we would argue that our analysis should be better 
suited to finding ongoing dust production by stars. On the 
basis of the above discussion, we advise caution when in- 
vestigating small infrared excesses in such cases and note 
the benefits of ph ase-matched, high-resolution infrared pho- 
tometry fsee alsolSloan et all l2010l ; iMcDonald et alj l2011bl ; 



iMomanv et alj|2012h 



5 CONCLUSIONS 

In this work, we have demonstrated the use of spectral en- 
ergy distribution fitting to determine the fundamental pa- 
rameters of the Hipparcos star sample. We have further used 
this information to quantify excess flux over the entire opti- 
cal, and near- and mid- infrared region of each SED. We have 
combined these excesses to determine those stars showing an 
excess of infrared flux, and cross-correlated literature iden- 
tifications to examine the cause of that excess over different 
regions of the H-R diagram, com paring our results t o the 
key studies of llta et all l|20ld ) and lGroenewegenl (|2012tl . We 
find we cannot reproduce the infrared excess and dust pro- 
duction claimed by the latter paper. 

Our analysis has focussed on the Hipparcos data cat- 
alogue: data which is now over 20 years old and, despite 
showing its age, provides the best estimate of distances to 
nearby stars we have. The launch of Gaia, and the comple- 
tion of further all-sky surveys such as Pan-STARRS, SDSS 
and WISE, will allow a similar analysis to be performed on 
many times more objects. Automated techniques, building 
on the kind demonstrated here, will be necessary to analyse 
and classify the objects which come from these surveys, in 
order to gain a full and comprehensive understanding of our 
corner of the Galaxy and its inhabitants. 
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APPENDIX A: REMOVAL OF BAD DATA 

We describe in this Appendix the sequentially-applied cuts 
we use to remove bad data from our catalogue. 



Al Cut #1 

The first cut was designed to remove extremely cool sources 
and sources with distinctly double-peaked SEDs from the 
sample. These tend to be heavily-enshrouded or heavily- 
extincted objects which we cannot accurately model. 

Stars were removed from the main catalogue if they 
had at least two mid-infrared (3.5-25 /im) bands brighter 
than all their optical/near-IR (w'-A" s -band) data. This cut 
removed 241 objects from the catalogue, of which 98 objects 
have no IJHK s -ba.nd data, leaving 109 436. 



A2 Cut #2 

The second cut was designed to remove unphysically low val- 
ues from the WISE data, which are much too faint to come 
from a stellar object detectable by Hipparcos. A WISE pho- 
tometric datum was removed from the combined catalogue 
if its flux was below 100 ^tJy. In this way, WISE photometry 
was deleted from 281 objects. 



A3 Cut #3 

The third cut acts to remove more underluminous W\ and 
Wi data. To assist in this, we define two fluxes, Fio and 
F20: Fio is the 10- (im flux defined by (in order of preference) 
the W 3 , AKARI [9] or IRAS [12] flux, and F 20 is defined 
similarly by the W 4 , AKARI [18] or IRAS [25] flux. 
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This cut applies to all stars where Fio > F20 and stars 
where no F\o measurement exists. This requirement pre- 
vents the selection of dusty sources where the SED reaches 
a minimum between 3 and 8 /im but where dust emission is 
insufficient to be picked up by Cut #1. W\ and/or W2 are re- 
moved if it their fluxes are less than 40% of both the 2MASS 
(or DENIS, where substituted) if s -band and the -Fio fluxes. 
Where no measure of -Fio is available, F20 is used instead. 

The cut removes WISE data from 273 objects. 

A4 Cut #4 

The fourth cut removes more underluminous Wi data. This 
cut detects objects where the flux in Wi is <40% of the flux 
in W2, provided the flux in W2 is: (1) non-zero, (2) greater 
than the flux in W3 and (3) less than the flux in K s (where 
W 3 and K s fluxes exist). The Cepheid X Sgr (HIP 87072) 
is excluded from this cut as a special case. This removes 38 
bad Wi datapoints from the catalogue. 

At this point, sources were also removed from the cata- 
logue if they had data in five or fewer of the observed bands. 
This removed 1479 objects, leaving 107 957. Most of the 
deletions were either detected by Hipparcos and IRAS, or 
Hipparcos and SDSS. 

A5 Cut #5 

The fifth cut is designed to remove single bad photometric 
data points. It operates on all filters shortward of 8.6 /im, 
except the Hipparcos / Tycho data. 

For each Hipparcos object, we determine the worst- 
fitting filter out of those listed in the previous paragraph, i.e. 
the filter with the greatest value of R or 1/R. If this value 
is five times greater than the next largest R or 1/R, it is 
removed. For example, if J and K s are the two worst-fitting 
filters, Rj = 51 and Rks = 10 then the J-band datum will 
be removed, whereas if Rj = 49 it will not. 

This cut removes 78 bad datapoints from the catalogue. 
Of these, 60 are DENIS /-band fluxes. 

A6 Cut #6 

The sixth cut was done manually to remove five stars where 
IR photometry is clearly confused due to blending, variabil- 
ity or background. These stars are HIP 60782, 80057 and 
88267 (all binaries), and HIP 82850 and 82611. This left 
107 952 unique catalogued objects. 

A7 Cut #7 

The seventh cut removes MSX B\ and/or B2 data when 
these data have a higher flux than than the K s filters and 
one of the W3, AKARI [9] or IRAS [12] filters. While in 
principle this could remove points from SEDs peaking be- 
tween 2.2 and 12 /mi, no objects seem to be affected by this. 
Dustless stars tend to peak at wavelengths shorter than 2.2 
/im, while dusty stars with double-peaked SEDs tend to have 
their second peak at wavelengths longer than 12 /im. Data 
was removed from 244 objects, though this affected only 
fraction of these, as MSX data is only used when WISE 
data is unavailable. There were 13 objects removed during 



this stage because they had insufficient photometry, leaving 
107 939 unique objects. 

A8 Cut #8 

The eighth cut is designed to remove underluminous SDSS 
data from saturated sources. SDSS g- and/or r-band data 
are removed if they are recorded to have less flux than the 
Hipparcos and Tycho Bt and Vr fluxes. This affects 7058 
objects, though most of these were already rejected in H2.ll 
An additional 224 points were manually removed from 
147 objects, mainly consisting of errant W\ and W2 fluxes. A 
further 39 objects were rejected for having insufficient data, 
leaving 107 900 unique objects. 

A9 Cut #9 

The ninth cut repeats cut #5, removing datapoints which 
have a goodness-of-fit more than a factor 2.5 (instead of 
5) worse than the next worst-fitting point. This cut was 
performed in two interations, the first removing a point from 
647 objects, the second removing a point from 59 objects. A 
final 120 points from 112 objects were edited out by hand. 
Following this, objects with less than five photometric points 
at <8 /im were rejected, leaving 107 619 unique objects in 
the final catalogue. 
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